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ABSTBAJT 



Becursive image segmentation with hierarchical scope 
views is a new technique in image processing which systemat- 
ically divides an image into smaller and smaller quadrants 
with each region within the quadrants having a structured 
descriptor. The regions are then processed to remove noise 
and to check for connected regions across the quadrant 
boundaries. The quadrants are then brought back together 
with the obscured target enhanced and distinguishable from 
the background. Infared image data of five different ship 
targets with the associated noise and inteference was 
processed by this technique with the target information 
being greatly enhanced. The procedures developed to eval- 
uate the data were found to be inadequate for the task, 
neccessitating hand evaluation to extract the target. The 
technique was proven to be a viable solution for extracting 
target information from infared data, but very slow 
processing times and inadequate evaluation procedures limit 
the usefulness of the program in its present form. 
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I. INTRODUC TIO N 



Many different approaches exist in the field of image 
processing to extract meaningful information from such 
sources as aerial photographs, infared images, or TV camera 
images. The major procedures of a typical image processing 
system are digitization, preprocessing, segmentation, 
feature extraction, and understanding. As presented in 
[Bef. 1], the purpose of segmentation is to partition the 
image into meaningful regions which depend on the problem 
being considered. In aerial reconnaissance, for example, a 
practical application of image segmentation may be to 
extract regions corresponding to terrain, industrial sites, 
or ships of interest. 

There are two ca tagories of approaches for the segmenta- 
tion process. The first category deals with methods which 
are based on examinirg an image on a point-by- point basis, 
called a point dependent process. An example of a point 
dependent process would be gray scale thresholding. The 
second category, referred to as a region dependent process, 
deals with techniques which utilize the image formation in a 
prescribed neighborhood. Examples of this category of 
approach are based on edge detection, boundary finding, or 
region growing which attempt to group points with similar 
characteristics into regions. 

A new method of segmentation utilizing methods of both 
catagories is by gray scale thresholding at hierarchical 
scope views [Ref. 2], which breaks an image into small quad- 
rants with the regions within the quadrants having a struc- 
tured descriptor and region number. This method uses many 
of the techniques first proposed by Ohlander (see [Bef. 3] ) 
and then formulated into the MOOSE program by Shafer 
[Ref. 3]. 
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Additional segmentation procedures are then applied to 
these small regions to pull out details and associate 
regions on oppsite sides of the boundaries. The segmented 
pieces are then processed to ascertain whether the region is 
valid or noise. All the pieces are then brought back 
together with the meaningful objects more prominently 
displaced. 

Previous segmentation processes have not provided good 
results when using infared data as the input. Interference 
caused by such variables as temperature differences, clouds, 
and target variations have caused the segmentation processes 
to fail in one way or another. A new program, called QUAD 
SPLIT, will be tested with infared data input of ships in an 
attempt to produce meaningful and useful information. To 
extract the information, new procedures will be required to 
take the output from the QUAD SPLIT program and process it 
into a form that can te readily understood. The formulation 
and testing of these new procedures, and the testing of the 
QUAD SPLIT program utilizing infared data of ships as the 
input is the basis for this research work. 

An overview of the recursive image segmentation methods 
utilized and their limitations are given in Chapter II, 
while a more detailed treatment of recursive image segmenta- 
tion at hierarchical scope views is presented in Chapter 
III. The problems associated with the program QUAD SPLIT 
and the development cf new procedures to remedy these prob- 
lems are the subjects of Chapter IV. The testing of the 
revised program with infared image data and the evaluation 
of the processed data is covered in the final chapter. 

The evaluation cf the processed infared data confirmed 
that recursive imace segmentation at hierarchical scope 
views is a viable technique for extracting target informa- 
tion from noisy infared images. The procedures developed 
for the evaluation proved to be inadequate which required 
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that visual and hand calculations be conducted to effec- 
tively extract the target. These inadequacies could be 
easily remedied since the required information is present in 
raw form and only requires proper processing. A major 
concern noted during testing is the processing time required 
for each image. For complicated images with numerous guad- 
rant segmentations and boundary checks, processing time can 
te in excess of twenty minutes. 



II. RECURSIVE IMAGE SEGMEN TAT ION METHODS 



This chapter is provided to give a brief overview cf two 
of the methods used in this work for the segmentation cf 
image data. The first section is an introduction to the 
MOOSE program and seme of its limitations. The second 
section then discusses a method by which two of the MOOSE 
limitations can be remedied. 

A. RECURSIVE SPLITTIBG 

1 . As sumpt i ons 

MOOSE is a program which implements a segmentation 
procedure using the assumptions proposed by Ohlander. The 
primary assumption for this technigue is that the surfaces 
in the image will be represented by connected sets of pixels 
with a very close gray scale measure. The data used 
throughout this work are infared images represented as a 0 
to 256 gray scale over a plane. The following additional 
assumptions are related to overcoming some of the flaws in 
the primary assumption: (1) Each of the connected sets of 
pixels, forming a patch, is assumed to produce a gaussian 
peak when histogrammed over the gray scale; (2) The histo- 
gram of a collection of these patches is assumed to clearly 
indicate a separation between peaks, with the relative 
minima of the histogram lying between these peaks rather 
than cutting peaks in the middle. 

2 . Alg orith m 

Ihe basic idea of the algorithm is to split the 
image into regions, then split these regions into smaller 
regions, and so on, until only very small regions remain. 
The structure is thus recursive. 
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one 



An image is considered as consisting of only 
region at the beginning. If the region area is larger than 
a constant minimum size criterion, it will be segmented 
further. A histogram of the region is then calculated and 
the shape of the histogram is analyzed to yield a set of 
peak intervals. The relative minima between the peaks is 
then used to produce good thresholds levels. The region is 
then divided by these threshold levels into candidate 
patches that are each assigned a distinct numeric score. 
The patches are then collected by a connecting region proce- 
dure to reveal their geometric relationship. A patch must 
have an area greater than some constant minimum noise area, 
otherwise it is considered to be noise and merged with the 
surrounding region. The patches are then stored in a region 
record. This is then repeated for other regions until no 
further regions require splitting. the procedure and asso- 
ciated criteria are listed as follows: 

Procedure Criteria 



Fetch a region Minimum area 

Histogrammi ng 

Peak selection 

Thresholding 

Connected patch finding 

Noise elimination Minimum area 

Region collection 
Repeat until region can not 
be split further. 

3 . limi tation s 

This recursive method has been applied yielding 
better results than the simple thresholding of gray scale 
approach. But there remain four main limitations: 



a. Majority Buie Problem 

The majority rule problem occurs when large and 
small objects (or areas) occur in the same region. Large 
objects will show up strongly in the histogram and will 
dominate. The smaller objects may not cause strong peaks in 
the histogram and thus may be covered by the large distribu- 
tional spread of the larger objects. Therefore, these 
smaller objects will not appear in the segmented output. 

t. No Well Eefined Peaks 

If there are no well defined peaks in the histo- 
gram then this procedure is ineffective. An image full of 
many small objects (i.e., cars in a parking lot) has a 
histogram which tends to have a relatively broad and flat 
distribution with no dominate peaks or valleys. Therefore 
these objects will be lost in the output. 

c. Region Combination 

Regions selected in early stages cannot be 
combined in later stages. The region boundary formed in the 
early stages will never change later on. Only more new 
boundaries will show up inside the old boundary. 

d. Constant Gradient 

Regions with constant intensity gradient cannot 
be detected. This can result in the splitting of an object 
that has a constant gray scale gradient instead of recog- 
nizing it as one object. 
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B. BECOBSIVE SPLITTIBG AT HIEBABCHICAL SCOPE VIEWS 



1 . Gen era l ities 

This approach was developed and implemented to avoid 
the first two disadvantages of MOOSE. The MOOSE algorithm 
is incorporated into the program with additional tests and 
procedures added. It uses a hierarchy of scope views with a 
large scope view (large area) at high levels and snail scope 
views (small areas) at lower levels. A quad tree (4 branch) 
development is used for splitting. Starting with a 256 X 
256 pixel representation, the image is divided into 4 guad- 
rant scope views by recursive image splitting. A test is 
then performed on the guadrant scope view to determine if 
further splitting is required. If the test results are 
acceptable, the quad tree development for this guadrant 
scope view stops at this point and is called a terminating 
node. If the test is not satisfactory for the. scope view, 
then it is further divided into four more quadrant scope 
views. Each of these is tested in turn which will either 
terminate or again he split into four guadrant scope views 
depending on the test result. This continues until a 
minimum area criteria is violated for splitting or all scope 
views test satisfact orily and terminate. 

The test used to check for acceptable results is 
based on the spread width of the histogram generated for the 
scope view. An object containing a large number of objects 
and a wide variation cf gray scale levels would tend to give 
a broad distribution histogram with no well defined peaks. 
If the spread width exceeds a maximum value set up by the 
test, then the test is considered not satisfactory and the 
scope view will be split into guadrant scope views. The 
test is then performed on these smaller scope views. Thus 
because the larger quadrants are broken up into smaller 
quadrants with a better possibility of finding well defined 
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histogram peaks, the majority rule and no well defined peaks 
problems associated with MOOSE are relaxed. 

2 . Bou ndar y Pr ob lem s 

Cne of the major problems of guad tree generation is 
that discontinuous boundaries may exist from one scope view 
to another. Two neighboring scope views with adjacent 
borders sometimes have different peak interval values which 
may lead to a non-clcsed region boundary in one scope view 
being discontinuous across the border. A boundary checking 
procedure was developed to alleviate this problem and is 
discussed in greater detail in the next chapter. 

C. CCHCLODIHG BEMABKS 

Two very related methods have now been introduced for 
the segmentation of images. It can be seen tha the recur- 
sive splitting at hierarchical scope views is just a refine- 
ment of the basic MOOSE program with the program QUAD SPLIT 
using MOOSE as its major building block. A more detailed 
discussion is provided in the following chapter. 
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III. BBCUBSIVE SPLITTING WITH HIE B ABCHICA L VIEWS 



In order to take infared data of a ship and process it 
into a useful output requires extensive preparation and 
knowledge of the system used. This chapter highlights the 
different problems involved in making these preparations and 
a more detailed look into the workings of the major proce- 
dures of the segmentation program used for this work. 

A. GEHEEAL 

The QUAD SPLIT program is the primary tool that was used 
in this research. The MOOSE algorithm is the basis fcr the 
QUAD SPLIT program. Therefore a thorough understanding of 
all the procedures of both was required. The version of the 
QUAD SPLIT utilized was written in PASCAL and contains 
approximately 5000 lines of code. Thus a more than tasic 
knowledge of PASCAL was also required to ascertain what was 
happening in each of the many procedures and where all the 
information was stored for future use. With no prior knowl- 
edge cf this language before the start of the research, very 
much valuable time was consumed in acquiring the needed 
knowledge as new complexities, structures, and formulations 
were encountered. 

QUAD SPLIT is a very complicated program which was 
implemented without regard to speed of execution. This made 
the tracing through and understanding of the inner workings 
of the program very difficult because many of the procedures 
are nested within other procedures which again are nested in 
still ether procedures. 

The infared data available for evaluation and testing is 
on magnetic tape in a 64 X 256 array of gray scale levels 
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from 0 to 255. QUAD SPLIT requires an input data file of a 
256 X 256 array. This data file was generated by formu- 
lating a short PASCAL program to copy the infared data in 
the first 64 lines cf the array and putting the value 1 in 
the rest of the array. With this type of input file 
supplied to QUAD SPLIT, it required an average of 10 minutes 
for the program to run to completion and produce an output 
file. 

B. MAJOE EBOCEDURES 

Hie rarchical S egmentatio n 

QUAD SPLIT takes the input data and reads it into a 
global record which is available throughout the program. 
The 256 X 256 array is considered the initial scope view and 
is at the highest level of the hierarchy at level 8 (2 ** 8 
= 256). When a scope view is divided, it drops one level 
(i.e. 128 X 128 => level 7, 64 X 64 => level 6, etc.). The 
level 8 image is automatically divided into four quadrant 
scope views at level 7. The upper left-hand quadrant, 
called NW, scope view (from 0,0 to 128,128) is tested first. 
The histogram for this scope view is calculated which is 
then checked for significant peaks. The candidate peaks are 
then tested for spread width to determine if further divi- 
sion is required. If so, the level 7 scope view is divided 
into four level 6 quadrant scope views. The NW scope view 
is again checked first to determine whether further division 
is required. A minimum area criteria for scope view area 
precludes division below level 5 (32 X 32) . 

When a scope view requires no further division, the 
local ninima between peaks is determined to set up threshold 
levels for assigning unique region numbers to the position 
in the scope view that fall within the threshold levels. 
For example, assume there are peaks in the histogram at gray 



18 



scale levels 100, 175, and 220with local minimas at 150 and 
200. Then for each array position in this scope view with a 
gray scale level below 150, a unique region number would be 
assigned (i.e. 4), while array positions for gray scale 
levels between 150 and 200 would be assigned another unique 
region number (i.e. 132), and for array positions with gray 
scale levels greater than 200 would be assigned yet another 
region number (i.e. 777). Thus, three distinct regions 
within the scope view have been created. The scope view has 
now satisfied its initial tests and is considered terminal. 

The next scope view to be evaluated is the upper 
right-hand, called NE, scope view on the same level as the 
last completed NW sccpe view. This continues on tc the 
lower left-hand (SW) scope view and finally the lower right- 
hand (SE) scope view. For example, start at level 8 which 
automatically goes tc level 7. The NW scope view of level 7 
is checked first and assume nc further division is reguired. 
The NE scope view of level 7 is then evluated, and assume 
this reguires division to level 6. Here the NW scope view 
is evaluated first and assume it is terminal, followed by 
the NE scope view (also terminal) , followed by the SW scope 
view (terminal) , and finally the SE scope view which is 
assumed to require further division to level 5. Assume all 
four scope views here at level 5 are terminal, which causes 
the SE scope view at level 6 to become terminal, which 
completes the NE scope view at level 7 and causes it to be 
terminal. This same sequence is applied to the other scope 
views at level 7 and at any lower levels that may be 
reguired until all the scope views at level 7 are terminal 
which in turn causes level 8 to be terminal. This example 
is depicted in Fig. 3.1. 
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2 . 



Noise El imin a tio n 

The level 8 array now has- a version which consists 
of the different regicn numbers. The number of times each 
region number occurs in the array gives the area for that 
specific region. This area is then compared to a constant 
which represents the minimum area a region must have in 
order not to be considered noise. If the region is less 
than this minimum area it is eliminated and is assigned the 
region number that surrounds it. The area for this 
surrounding region is then updated and the check continues 
until all regions have been checked. 

3 . Boundar y Che c kin g 

The array of region numbers could now be made avail- 
able in output to give a segmented representation of the 
original infared data. Problems may occur at the boundaries 
of the different scope views because the histogram and 
threshold levels used may be different in the adjacent scope 
views. Take for example the simple situations presented in 
Pig. 3.2. In scope view A, histograming and thresholding 
produce three regions that end at the boundary. In the 
adjacent scope view, scope view B, histogramming and thresh- 
olding may produce only two regions at the same boundary 
which may or may not coincide with the regions of the first 
scope view. 

In scope view D there are two regions in which one 
region is partially enclosed by the other with an abrupt 
ending at the scope view boundary. On the adjacent 
boundary, scope view C, there are no breaks in the boundary 
due to different results for histogramming and threshold 
levels, even though it appears that there should be a 
continuation of the partially enclosed region of scope view 
D across the boundary. 
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A procedure called SPLIT CHECK was thus incorporated 
into QDAE SPLIT to eliminate some of these potential prob- 
lems. The scope views into which the image was hierarchi- 
cally split is recalled. The four guadrant scope views at 
the lower levels are checked first. The rightmost column of 
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the NW and SS scope views, called the NS boundary because it 
is vertical and runs from top to bottom or north to south, 
is used to start the check. Each position of this column 
has its region number compared to the region number of the 
position just below it. If the two values are different 
then the boundary between two regions within the scope view 
has been detected and is called a break. When a break is 
detected, the row cf the position where this occurred is 
remembered. The adjacent scope view boundary (leftmost 
column of the NE and SE scope views) is then checked within 
a range cf +5 to -5 positions of the row in which the break 
was found. The check performed is the same as detailed 
above for finding a break. If a break is found within this 
range, then there is no problem for the desired result to be 
realized. If no break is found on this adjacent boundary, 
then the adjacent scope view is divided again, if possible, 
with histogramming, thresholding, and region number assign- 
ment again taking place. Boundary checking is then recom- 
menced with the hope of now finding a break in the 
appropriate area on this adjacent boundary. 

If this still fails, then the thresholds of the 
original scope view regions for which the break was found 
are recalled and averaged. If this new value falls within 
the peak interval for thresholding of the adjacent scope 
view region, then it is possible to dissect the peak 
interval into two new intervals and thus create two regions 
and a possible break. This is called threshold propagation 
and is depicted in Fig. 3.3 . If this still fails to 
produce a break, then nothing more is attempted and the 
boundary checking procedure continues. 

This checking for a break is done at all positions 
down the column except at rows divisible by 32 where scope 
view boundaries exist and thus a break would occur. After 
this column is checked, its adjacent column is likewise 
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Boundary 




Region C 



No Neighbor Break E xists 

Assume that the peak interval of Region A >= the peak 
interval of Region B. 

New threshold for region C = 1/2 (min. of interval A 
max. of interval B) . 

Cl > new threshold C2 < new threshold 




Ne ighbo r Break N ow E xists 

. — . — 



Figure 3.3 Threshold Propagation. 
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checked and compared with the original boundary column if a 
break is detected- Next the bottom most row of the Nfl and 
NE scope views, called the EW boundary is used to check for 
breaks in the column positions. Any breaks found are 
checked in likewise fashion across the boundary. This is 
again done for the ether side of the boundary until all 
boundary positions inside the four scope views have been 
checked. This same procedure is repeated for all levels 
until the four scope views at level 7 have been finally 
checked, and the procedure ends. 

4 . Cut put G ener a tio n 

Region numbers have now been assigned to all the 
regions in the different scope views and are available in a 
256 X 256 array. To get a hard copy of this region map from 
a printer, the region numbers must first be converted to 
single letters, numbers, or symbols so that each array posi- 
tion representing the different regions can be represented 
by a single character- This is accomplished by using 55 
different characters from the keyboard. The region numbers 
are first put through a modulo 55 function and the remainder 
is correlated to one of the 55 characters. This array of 
characters is then stored and is available for output to a 
printer in four 128 X 128 pieces. 

By linking the QUAD SPLIT program with the PLOT 10 
program available on the VAX-750, an output depicting the 
edges of all the regions can be displayed on the TEXTRCNICS 
console. 
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IV. IMAGE RESTRUCT URIN G 



A. LIMITATIONS OF QOAD SPLIT 

1 . Eeqion Numb e ring 

As detailed in the previous chapter, each region in 
the scope view has its own region number which correlates to 
its own region character in the hardcopy output. At bound- 
aries there will be different characters, even though both 
sides of the boundary may have the same gray scale level in 
the original infared data. Therefore background features 
such as sky, horizcr, water, etc., are broken up and 
displayed by many characters depending on how many times the 
image has been split ever that area. This problem is some- 
what eliminated by useing the PLOT 10 program in some ether 
form cf output because only the edges are displayed. lhe 
output must then be visually interpreted to find the desired 
features . 

2. Eeq ion Map Ou tpu t 

Problems of duplicate characters can arise in the 
output because of using the modulo 55 function to operate on 
the region numbers when assigning characters to the listing 
output. In QUAD SPIIT the maximum number of regions that 
can he generated is limited to 64, while region numbers can 
have values from 1 to 4096. Thus duplication of characters 
is highly likely, but their occurence in adjacent regions or 
across boundaries is not often. 

3 . Back groun d Fe atures 

The procedures and tests incorporated into QUAD 
SPLIT for the calculation of the area and the centroid of 
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each individual region, and for making a selection of good 
features are limited because of the region propagation 
problem. Background features that may occur in the image 
could have been broken up at least at scope view boundaries. 
Thus, area and centroid information become inaccurate for 
selection of a target from the background regions. 
Additional procedures for more extensive geometry comparison 
is reguired to distinguish a target from the rest of the 
regions. 

B. AIGOBITHH FO BHULATIOH 

1 . Beg ion Number Prop ag at ion 

"Various ideas were considered to accomplish the 
aesiered region number propagation across the boundaries. 
Cne method considered was to count the number of times a 
region number occured at a boundary and compare the result 
with a similar count on the opposite side of the boundary. 
The assignment of the same region number to those that 
occured almost an equal number of times on both sides of the 
boundary would then take place. This was quickly rejected 
because the count on one side of the boundary may occur at 
the top while on the opposite side of the boundary it could 
occur at the bottom, and thus the same region number could 
le assigned to completely unrelated areas. 

Another idea- was to count the number of times two 
pairs of region numbers occured across a boundary and to 
assign the same region number to the pairs that occured most 
frequently or that exceeded a set minimum count. This was 
also rejected because one side of a boundary may have only 
one cr two region numbers while the opposite side of the 
boundary could have numerous region numbers and when counted 
in pairs could all realize the same result. Then which 
region numbers should propagate and which should not would 
have to he decided. 
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Ihe idea that was finally selected as the most 
likely to produce the desired result was to determine where 
breaks were encountered on both sides of the boundary and 
then assign the same region numbers to the pair of regions 
that occur across the boundary on either side of the break. 
If a break is detected on one side of the boundary but with 
no adjacent boundary break found within a specified range, 
then all regions would retain their present region numbers. 
This is the basis upon which the region number propagation 
algorithm will be implemented. 

2. New Reg ion Map 

Since there are a maximum of 64 regions available 
from the QUAD SPLIT program, this implies that there are at 
most only 64 different region numbers utilized. Thus, by 
adding nine more characters to the list of 55 characters now 
available, all the 64 different regions could be assigned a 
unique character corresponding to its region number without 
duplication being a problem. 

3. Feat ure Tes ts 

Easic FORTRAN procedures for calculating area, 
perimeter, center of gravity (or centroid) , size, and 
compactness were known to be available on the VAX-750 system 
by using the SPIDER programs [Ref. 4J. These procedures 
would use the information generated from the region number 
propagation procedure to produce their respective results 
for each region. The retrieval of the desired targets could 
then be based on this information. 
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C. ALGORITHM I HPLEM E STATION 



1. B ackground R egion Chec k 

The implementation of the algorithm for region 
number propagation was based on the boundary check procedure 
of the QUAD SPLIT program. Since it was to help detect the 
background portions of the image, the procedure was named 
Background Region Check , or 3AKREG_CK for short. The 
procedure was to read the region number information from 
QUAD SPLIT into a new array which then can be changed 
without affecting the original record. 

The same type of boundary checking as discussed 
previously was utilized to detect a break. If a break was 
detected on the opposite side of the boundary within +5 to 
-5 positions of the original break, then the pair of two 
adjacent regions above (or to the left) of the break would 
be assigned the region number of the region to the left for 
NS boundaries (or the region on top for EW boundaries) , and 
the pair of adjacent regions on the other side of the break 
would be assigned the region number of the region below and 
to the left of the break for NS boundaries (or above and to 
the right of the break for EW boundaries) . A pictorial 
diagram of the above situation is given as Fig. 4.1. If no 
break is found in the adjacent boundary, then no changes are 
made. 

Cne potential problem that was realized at this 
point was the situation where multiple breaks on one side of 
the boundary were matched by offset multiple breaks on the 
adjacent boundary within the +5 to -5 position range. For 
example, (see Fig. 4.2), assume we are checking column 64 
and a break is encountered at row 15, with the region above 
the break designated region ’A' and the region below the 
break region •C’. A corresponding break on the opposite 
side of the boundary (column 65) occurs at row 19, which is 



29 



*1 



MS B oundary „ EH B oundary 





BE FOR E 





AFTE R 



Figure 4. 1 Example of Boundary Checking. 
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within the +5 to -5 position range. Call the region above 
this break region *B* and the region below as region *D*. 
Since corresponding breaks on both sides of the boundary 
have been found, then region propagation occurs with regions 
•A* and 1 B ’ being assigned region number ’A', and regions 
*C* and 'D* being assigned region number 'C'. Checking now 
continues down column 64 for another break and assume for 
this example that it occurs at row 23 with region *C* above 
the break and region *E' below. Checking adjacent column 65 
would shew a break at row 19 again since it is still within 
+5 to -5 positions of this new break in column 64, while an 
additional break in column 65 at row 24 would be missed. 

lo alleviate this problem additional steps were 
incorporated so that once two adjacent breaks have teen used 
for region number propagation, then neither can be used 
again by additional breaks. Also it was noted that only one 
pass along a given set of adjacent boundaries is sufficient 
for this check because once one boundary is checked and 
compared to its adjacent boundary, then all of the possible 
connected regions will have been recognized and nc further 
information would he gained by checking the adjacent 
boundary against the original boundary. 

2. Begion Map Output 

Ihe assignment of new region numbers from 1 to 64 
was implemented by a reiterative routine which would run 
from 1 to 4096 checking the region number array at each 
position. When an original region number was found, it was 
reassigned a new region number starting with 1 and going up 
to the maximum number of regions utilized. Each time a 
position was reassigned a new region number a counter would 
be incremented and this count was compared to the maximum 
number of postions available (65536) and would stop the 
iteration once this was reached. 
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Figure 4.2 Boundary Checking Example for Multiple Breaks. 
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To simplify processing throughout this set of 
procedures, this reassignment of region numbers was incorpo- 
rated with the inputing of the information from the QUAD 
SPLIT region number record. 

An additional nine characters was also added to the 
list cf characters utilized to represent the region numbers. 
The characters were then assigned to respective region 
numbers after the region number propagation procedure was 
completed to give a hardcopy output available in four 128 X 
128 segments. 

3 . Ge o metry 

The procedures for area, perimeter, and center of 
gravity (centroid) were the only SPIDER procedures utilized. 
The procedures for size and compactness both called the area 
and perimeter procedures and then applied a simple equation 
to produce their respective outputs. Since the area and 
perimeter outputs were already available, the following 
equations were used to directly calculate the size and 
compactness : 

size = (2 X area) / perimeter 
compactness = (4 X PI X area) / perimeter 

D. CCHCLUDIHG REMARKS 

The program has now been thoroughly reviewed and new 
procedures have been formulated and implemented. The proce- 
dures as written are provided as Appendices A thru E. The 
whole package must new be brought together so that test runs 
with real data can be observed. The testing and evaluation 
of the program is the subject of the next chapter. 
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V. TEST AND E?ALO ATION 



The basic ideas behind all of the programs have been 
looked over extensively and the algorithm changes to these 
programs are implemented. The testing of all procedures 
with the infared image data is the subject of this chapter. 
The evaluation of the test data and conclusions reached are 
in the final sections. 

A. EXPEBIHEHTA1 HES01TS 

Testing was an on going effort during the entire 
research period. As each small program or procedure was 
written, it would reguire testing to verify that it gave the 
desired results. This section will only cover the testing 
of the major programs and procedures that, have been 
discussed in detail in prior chapters. 

1 . £0 AD SPLIT 

The QOAD SPLIT program had been extensively tested 
prior to this work and the limitations that were encountered 
are documented in Chapter IV. The average running time of 
the program for the infared image data was approximately 10 
minutes. Fig. 5. 1 is an example of a portion of the 
segmented output listing. This example makes it obvious 
that the output is very disjointed at the scope view 
boundary occuring at column 126. It would be very difficult 
to gather any useful information from the output in this 
form . 
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Figure 5-1 An Example of QUAD SPLIT Region Map Output- 
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Bac kgro und Region Chec k 

The new procedures added to the QUAD SPLIT program 
are the real essence for the testing conducted in this work. 
As a new procedure was added to the original program, it was 
tested to check if any further refinements would be 
required. 

a. Region Number Propagation 

The first portion of the background region 
checking procedure to be implemented was the region propaga- 
tion algorithm and the associated new region map output 
listing. The initial tests of these procedures led to the 
detection of the problem where the same region symbol was 
being utilized mere than once in the same region map. After 
the incorporation of additional steps to rectify this 
problem, a series of tests was performed using five 
different infared images as the input data. The total time 
for each program to run was approximately 16 minutes, with 
the fastest completed in 10 minutes and the slowest in 21 
minutes. The time required for each run was determined by 
how many times the image data was split and by how many 
regions were formed. The data with larger targets had less 
regions and therefore ran the fastest. The same portion of 
segmented output listing as presented in Fig. 5.1 but with 
region number propagation in effect, is given in Fig. 5. 2. 
As can be seen, the scope view boundary at column 128 has 
been virtually eliminated and background features, such as 
regions 'B, a, and J f , which propagate across the entire 
output section are easy to pick out. For each of the five 
different infared images used, successful region number 
propagation was observed. The background region map outputs 
are presented as Fig. 5. 3, Fig. 5. 4, Fig. 5. 5, Fig. 5. 6, 
Fig. 5.7, Fig. 5. 8, Fig. 5. 9, Fig. 5. 10, Fig. 5. 11, and Fig. 5. 12. 
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Figure 5.2 Region flap Output with Region Number Propagation. 
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Figure 5.4 Image A EacJcground Region Hap — col. 129 to 256. 
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Figure 5.6 Iaage B Background Region Map — col. 129 to 256. 
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Inage C Background Region Map — col. 1 to 128 
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Figure 5.8 Iaage C Eackground Region Map — col. 129 to 256. 



43 






COLUMN 




os o s 



Figure 5-9 Inage D Background Region Hap — col. 1 to 128. 
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figure 5.10 Image D Background Begion Map — col. 129 to 256. 
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b. Geometry 



The procedure for obtaining the area, perimeter, 
centroid, size, and compactness for each of the regions was 
tested next. Each of the individual FORTRAN procedures were 
tested separately by short routines with dummy input data to 
ensure that the desired result would be obtained and to 
verify that all required parameters were supplied correctly. 
Each of these were tested satisfactorily and then the area, 
perimeter, and centroid procedures were incorporated into 
the basic program. Since size and compactness can be 
derived from the area and perimeter procedures, it was 
decided to save some calculating time by using the eguaticns 
to generate the desired result. A sample listing of the 
data generated for the regions shown in Fig. 5. 2 is given in 
Table 1. All five infared images were processed with these 
added procedures and less than one minute of processing time 
was added to the overall run time respectively. 

B. E UA1DATI0H OF DATA 

All of the tests have now been completed. The data must 
now be evaluated to determine if the desired results of 
picking cut the target from the background and noise are 
achieved. Time limitations precluded the generation of 
procedures to incorporate the program for this purpose. The 
following evaluation was conducted by visual examination of 
the data and some manual calculations. 

1 . Ba c kground Re gion s 

After reviewing the background region map outputs, 
it could he readily observed that the background sections 
were clearly recognizable. They were regions that went 
generally across many boundaries and contained a major 
portion of the area. From the Geometry data files, it was 
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1 












1 








TABLE 


1 










Sample Geometry 


Data 




REGION 


AREA 


PERI- 


CENTROID 


SIZE 


COMPACT 






TIETYE 


555 — 551 




“5555“ 


1 


26 


18 


8 


97 


2.89 


1.008 


8 


48 


20 


35 


148 


4.80 


0.513 


9 


70 


46 


33 


138 


3.04 


0.416 


E 


3275 


348 


12 


167 


7.72 


0.057 


F 


25 


42 


3 


152 


1.19 


0. 178 


G 


15 


24 


9 


164 


1.25 


0.327 


E 


1981 


752 


14 


95 


5.27 


6.044 


I 


9 


16 


46 


156 


1.12 


0.418 


J 


1198 


22 


31 


160 


109 


31.100 


I 


34 


3 2 


43 


172 


2. 13 


0.417 


M 


123 4 


334 


23 


94 


7.39 


0. 139 


C 


1230 


356 


37 


65 


6 .31 


0. 122 


D 


3815 


16 


48 


164 


477 


187.30 


E 


190 


1 14 


4 


182 


3.33 


0. 184 


X 


35 


42 


31 


131 


1.67 


0.249 


noted that the 


regions with 


the 


largest 


areas and size 


correlated 


with 


the background 


regions noted 


visually in all 


five cases 


. Therefore, 


all of 


these regions 


could be easily 


eliminated 


in a 


search 


for the 


good 


target . 


In the example 


of sample 


data < 


given in 


Fig. 5. 2 


and 


Table 1 


the regions with 


symbols *B 


, H, 


J, H, 


0, and 


U* 


would then be eliminated. 



For the five infared images evaluated, this type of hack- 
ground region determination would eliminate one third cf the 
regions and approximately ninty percent of the area. 
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2 . 



Size vs Co mpa ctness 

A correlation between size and compactness was then 
attempted for the remaining regions. From all the data 
obtained, no useful correlation could be obtained. Even 
after looking at the background region map outputs and 
selecting likely candidates to check for any correlation, 
none could be found that would cover the target and get rid 
of the other regions. Using the output map of Fig. 5. 2 as an 
illustration, the regions with symbols '8, 9, and X* seem to 
be the most likely candidates for being the target repre- 
senting the ship. Checking the size and compactness data 
for these three regions from Table 1 yields a size range of 
1.67 to 4. 80 and compactness range of 0.249 to 0.513. Two 
other regions fall within this size range, three more within 
the compactness range, and one region falls within both 
ranges. Similar results were observed for the other four 
images. Mere information is therefore required to make a 
good decision. 

3. Di mens ions of Region s 

The dimensions of the regions were not incorporated 
in any of the procedures for the program under test. After 
reviewing all the background region maps, it was noticed 
that there were numerous long and thin regions. These 
regions were noted tc be usually close to boundary areas and 
are generally caused by some noise or distortion in the 
infared image. If these regions were eliminated, an addi- 
tional twenty five percent of the total regions could be 
removed . 

The dimensions for the regions remaining after back- 
ground region elimination for the example sample were 
obtained from the region map output and the results are 
tabulated in Table 2 . By reviewing all the data from the 
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five infared images, it was noted that if any dimension was 
of length 3 or less, or if the column lenth divided by row 
length did not fall within a range of 0.1666 to 6, then 
these regions could he assumed to be noise and eliminated. 
Therefore, from Table 2 data, region symbols 'G, I, and L' 
could be eliminated. 



TABLE 2 

Sample Region Dimensions 



REGION ROW NU MBER , . COLUMN NUMBER 



31FBCI 


MIN' 


— rTTNGTH 


MTN TTAX" 


T'ETTGT H 


1 


3 


7 


5 


125 


128 


4 


8 


33 


37 


5 


140 


156 


17 


o 

> 


31 


36 


5 


130 


148 


19 


E 


1 


5 


5 


146 


158 


13 


G 


9 


10 


2 


159 


170 


12 


I 


46 


4 7 


2 


154 


159 


6 


1 


42 


45 


4 


159 


183 


25 


ti 


1 


9 


9 


162 


204 


42 


X 


30 


34 


5 


125 


139 


15 



Another use cf the dimensions of the regions could 
be to use them in conjunction with the centroid data to 
locate regions which are in close proximity to each other. 
The infared data of a ship usually breaks the ship up into 
different regions because of the difference in intensity cf 
the infared radiation given off by different areas of the 
ship. These regions would all be close together and this 
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fact could aid in the selection of the good targets. a 
simple test was devised for taking the regions that are 
remaining and checking to find out if both the row and 
column values of their respective centroids are within a 
given fixed value. For the five infared images tested, a 
value of ten seemed tc work satisfactorily. For the sample 
data cf Table 1 this test would be satisfied for region 
symbols *8 and 9' and *9 and X’. Therefore, a good choice 
for the good regions representing the ship would be the 
regions designated by symbols *8, 9, and X 1 , which were the 
same regions selected for being possibly good targets from 
the background region map of Fig. 5. 2. 

This test was applied to the other four test images 
with reasonable results being obtained. A potential problem 
could be with an image of many small objects or ships, such 
as ships in a convoy, which would produce regions not close 
together resulting in the failure of this approach. 

C. CCHCIUSIONS 

The research work reported in this thesis had as its 
major objective the extraction of usable target information 
from infared image data using the recursive image segmenta- 
tion with hierarchical scope view technigue. The five 
infared images used in the testing were processed suffi- 
ciently to yield usable information. 

The new procedures of this chapter incorporated into the 
main program proved tc be inadequate for using its output to 
make a decision as tc which region is the target. By visual 
and manual processing this decision was shown tc be 
possible. Further research will be reguired to resolve 
these inadequacies. 

There are still many problems that must be solved before 
this type of procedure can be useful in the field. The 
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major problem is the processing time. An image that is 
small with considerable noise present could result in many 
regions being formed. This could cause the processing time 
to easily surpass twenty minutes. This is clearly unsatis- 
factory in any real time situation. A computer engineering 
research effort to reduce the complexities and excessive 
nests and loops could be of great benefit. 

The region number assignment problem that was dealt with 
in the new procedures of Chapter IV was not completely 
remedied. Regions that are unconnected but with the same 
region symbol assignment occured in all of the tests. The 
assigning of unigue region numbers must be approached from 
within the main program if this problem is to be resolved. 
This could possibly he addressed in conjunction with some 
other research effort utilizing the QUAD SPLIT program. 

An initial objective of the research was found to be 
overly ambitious once the complexities of the programs began 
to surface. The objective of producing a possible hardware 
realization of the procedure was therefore not pursued. 
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APPENDIX A 



EEOCEDOBE EAKBEGHAP 



PROCEDURE EAKREGMAP 

(* This procedure takes the region number 
output of QUAE SPLIT and assigns new 
region numbers from 1 to 64 in a new 
array. Region number propagation is 
performed, the new region map is output, 
and geometry calculations are made. *) 

TYPE 

EACKREGIONMAP = RECORD 

BRMAP_PIX : ARRAY {GROW_NC, GCOL_NO} OF GREG_NO 

END; 

VAR 

BEEGMAP : °BACKR EGIONMAP ; 

I, J, X, Y, Z : INTEGER ; 

END A : BOOLEAN; 

FOUND : ECOLEAN; 

BEGIN 

NEW (EEEGMAP) ; 

WITH BREGMAP 0 DO 

FOE I := 1 TO GEOW_MAX DC BEGIN 
FOR J := 1 TO GCOL_MAX DO BEGIN 

BRMAP_PIX {I, J} := GEEGMAP°.GRNAP_PIX {1 ,0} ; 
END; 

END; 

Y := 0; 

ENDA := FALSE; 

Z := 1 ; 
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FOF X := 1 TO 4096 DO BEGIN 
IF ENDA = FALSE THEN BEGIN 
FOR I := 1 TO 256 DC BEGIN 
FOR J := 1 TO 256 DO BEGIN 
IF BREGHAP°.BRMAF_PIX £I,J} 
IF Z > 63 THEN Z := 64; 
BREG MAP®. BRMAP_PIX{I,J} 
FOUND := TRUE; 

Y := Y ♦ 1; 

IF Y >= 65536 THEN ENDA 
END; 

END; 

END; 

IF FOUND = TRUE THEN BEGIN 
Z := Z + 1; 

FOUND := FALSE; 

END; 

END ; 

END; 

B A K B N D_C K (1) ; 

OUTBAKMAP; 

GECMETRY ; 

END; 



X THEN BEGIN 

z; 



TRUE; 
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APPENDIX B 
PEOCEDURE EAKBND_CK 

FROCEEURE EAKBND_CK (QDADNO : QNODE_NO) ; 

(* Background boundary check *) 

(* This procedure segments the image into 
guadrants and then calls the function 
BAKSPLT_CK to check the quadrant 
boundaries for breaks. *) 

VAE 

EDR_CLEAN, NS_BDE, FWAR D_CK : BOOLEAN; 
BER_SIZE, BDR_HA1E : 32. .256; 

EE ; GRON_NO; 

BC ; GCCL_NO ; 

BLEVE1 ; LEVEL_NC; 

3EGIN 

ELEVEI := QN {QUAENO} °. QLEVEL; 

WEI1ELN (LISTING, * QUAD NODE 1 , QUADNO) ; 

IF CN {QUADNO} °. STATUS = SPLIT THEN BEGIN 

IF CN {QN {QUAENO} °.NW_SCN} °. STATUS = SPLIT 
EAKBND_CK (QN {QUADNO} °. NW_SON) ; 

IF CN {QN {QUADNO} °.NE_SON} °. STATUS = SPLIT 
EAKBND_CK (QN {QUADNO} °. NE_SON) ; 

IF CN {QN {QUAENO} °.SW_SCN} °. STATUS = SPLIT 
EAKBND_CK (QN {QUADNO} °. SW_ SON) ; 

IF CN {QN {QUAENO} °.S E_SCN} °. STATUS = SPLIT 
EAKBND_CK (QN {QUADNO} °. SE_SON) ; 

END; 

NS_EEB := TRUE; 

FWARE_CK := TRDE; 

BE := CN {QUADNO} °. GROW ; 



THEN 

THEN 

THEN 

THEN 
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BC := QN [QOADNO} 0 . GCOL ; 

BDR_SIZE := 2 ** EIEVEL; 

BDR_HALF := BDR_SIZE DIV 2; 

EER_CLEAN := FALSE; 

R BILE BDR_CLEAN = FALSE DO BEGIN 
EDE_CLEAN := T ROE ; 

WRITELN (LISTING, 'NS FORWARD CHECK*); 

IF B AKSPLT_CK (NS_BDR, FRARD_CK) THEN 
EDR_CLE AN := FALSE; 

ENE; 

WRITELN (LISTING, * NS REVERSE CHECK'); 

IF BAKSPLT_CK (NS_BDR, NOT FWARD_CK) THEN 
EDR_CL EAN := FALSE; 

END; 

RRITELN (LISTING, *ER FORWARD CHECK') ; 

IF B AKSPLT_CK (NOT NS_BDR , FH ARD_CK) THEN 
EDR_CLE AN := FALSE; 

ENE; 



WRITELN (LISTING, ’ EW REVERSE CHECK'); 
IF BAKSPLT_CK (NOT NS_BDR ,NOT FRARD_C K) 
EDR_CLEAN := FALSE; 

ENE; 

IF BDR_CLEAN = EALSE THEN BEGIN 

IF QN (QN {QOADNO} 0 . NW_SON} °. STATUS 
3 AKBND_CK (QN {QOADNO} °. NR_SON) ; 

IF QN {QN {QUADNO} °. NE_SON} STATUS 
BAKBND_CK (QN {QOADNO} °. NE_SON) ; 

IF QN {QN {QOADNO} 0 . SR_SON} °. STATUS 
BAKBND_CK (QN {QOADNO} °. SR_SON) ; 

IF QN {QN {QOADNO} 0 . SE_SON} °. STATUS 
BAKBND_CK (QN {QOADNO} °. SE_SON) ; 
END ; 

ENE; 

END ; 



THEN 



= SPLIT THEN 
= SPLIT THEN 
= SPLIT THEN 
= SPLIT THEN 
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APPENDIX G 



FUNCTION BAKSPLT_CK 

FUNCTION B AKSPLT_CK (NS : BOOLEAN; 

EORWARD : BOOLEAN) : BOOLEAN; 

(* This function checks each boundary within 
a quadrant for breaks on both sides of the 
boundary. If this occurs, the adjacent 
regions across the boundary are assigned 
the same region number. This is referred 
to as region number propagation. *) 

TYPE 

BBK_TYEE = 1.. 257; 

VAR 

PEEV_TEE, PBEV_EEG, NOW_REG, NNOW_REG : REG_NO ; 
SC, NSC : GCOL_NO ; 

BC_LIMIT, BRK, BB_LIMIT : ERK_TYPE; 

SB, NSB : GBOW_NC; 

NERK, PREV_N BEK : -3. .255; 

ATEE : BOOLEAN; 

I, J, K : INTEGER; 

BEGIN 

E AKSPLT_CK := FALSE; 

E RE V_NERK : = 0; 

PEEV_TEE := 0; 

WITH EREGHAP 0 EO EEGIN 
IF NS THEN BEGIN 

IF FORWARD TEEN BEGIN 

SC := BC + EDR_HALF - 1; 

NSC := BC+ EDR_HALF ; 

END; 
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ELSE 3EGIN 

SC := BC + EDR_HALF; 

NSC := BC+ EDR_HALF - 1;' 

END; 

SR ;= BR; 

END; 

ELSE BEGIN 

IF FORWARD TEEN BEGIN 

SR ;= BR + EDE_HALF - 1; 

NSR ;= BR + EDR_HALF ; 

END; 

ELSE BEGIN 

SR := BR + EDR_HALF; 

NSR := BR + EDR_HALF - 1; 

END; 

SC ;= BC; 

END; 

PEEV_REG ;= BRMAP_PIX (SR, SC}; 

NCW_REG ;= PREV_REG; 

IF NS THEN BRK := BR ELSE BRK := BC; 

BB_IIHIT : = BR + EDR_SI2E; 

BC_LIM IT := BC ♦ EDR_SIZE; 

WHILE {(NS AND {EEK < BR_LIMIT) ) OR ((NOT NS) 

AND (BRK < BC_LIMIT) ) ) 

EO BEGIN 

WHILE (({NS AN! (BRK < BR_LIMIT) ) OR ((NOT NS) AND 
(BRK < BC_LIMIT))) AND ( NOW_REG = PREV_REG) ) 
LC BEGIN 

PREV_REG :+ NCW_REG; 

ERK := BRK + 1; 

IF (BRK < 257 ) THEN IF NS THEN 
N 0 W_ R E G := ERHAP_PIX (BRK, SC] 

ELSE NOW_REG := BRMAP_?IX (SR, BRK} ; 

END; 

IF {((BRK - 1) MOD 128) <> 0) 
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AND (((BEK - 1) MOD 64) <> 0) 
AND ( ( (BR K - 1) MOD 32 <> 0) ) 
THEN BEGIN 



ATEE := TROE ; 

NBRK := BR K -4; 

IF NS AND (NBRK < BR THEN NBRK := BR ; 

IF (NOT NS) AND (NBRK < BC) THEN NBRK := BC; 

FOR I := 1 TO 8 DO BEGIN 
IF NS THEN BEGIN 

NP REV_REG := BRMAP_PIX (NBRK, NSC}; 

NNOW_REG := BR MAP_PIX (NBRK + 1, NSC}; 

END; 

ELSE BEGIN 

NPREV_REG := BRMAP_PIX (NSR, NBRK} ; 

NNOW_REG := BRMAP_PIX [NSB, NBRK + 1}; 

END; 

IF ( NPREV_REG <> NNOW_REG) THEN BEGIN 
ATEE := FALSE; 

IF FORWARD THEN BEGIN 

IF (PEEV_NBRK < NBRK) THEN BEGIN 
IF NS THEN BEGIN 

FCR K : = BR TO (BR_LIMIT - 1) DO BEGIN 
FOR J : = NSC TO BC_LI MIT-1) DO BEGIN 
IF BRM AP_PIX (K , J} = NPREV_REG 
THEN BEGIN 

BRMAP_PIX (K, J} ;= PRE V_REG ; 
END; 

IF BRMAP_PIX (K, J}= NNOW_REG 
THEN BEGIN 

BRMAF_PIX (K, J} := NOW_REG; 

END; 

END; 

END; 

END; 

ELSE BEGIN 
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FCR K := NSR TO (BR_LIMIT-1) DO BEGIN 
FOR J := SC TO (BC_LIMIT-1) DC EEGIN 
IF BRMAP_PIX (K, J} = ?REV_REG 
THEN BEGIN 

BR MAP_PI X {K , J} := PREV_REG ; 
END; 

IF BRM AP_PI X {K,J} = NNOH_REG 
THEN BEGIN 

BRMAP_PIX (K, J} := NOH_REG; 

END; 

END; 

END; 

END; 

PRE\Z_NBRK := NBRK ; 

END; 

END; 

END; 

ELSE BEGIN 

NBRK ;= NBRK + 1 ; 

IF NS AND (NBRK < 3R) THEN NBRK ;= BC; 

IF (NOT NS) AND (NBRK < BC) THEN NBRK: = BC ; 
IF NS AND (NBRK > BR_LIMIT - 2) THEN 
NBRK := BR_LIMIT - 2; 

END; 

IF (NOT NS) AND (NBRK > BCJLIMIT - 2) THEN 
NBRK := BC_LIM IT - 2; 

END; 

END; 

END; 

IE ATEE THEN BEGIN 

WRITELN (LISTING,’ EREAK AT =',BRK, 

’NEIGHBOR AT’ , N3RK) ; 

PREV_TEE := BRK; 

END; 

END; 
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f IEV_SEG 
ENE; 

END; 

ENE ; 



NCW_REG; 











APPENDIX D 
EBOCEDOBE GOTBAKMAP 

PROCEDURE OUTBAKMAP (* Background Region flap Output *) 

(* This procedure assigns symbols to the 
region numbers for printer output. The 
output is formatted into four pieces of 
128 X 128 tc produce the entire output. *) 

VAR 

ECH ; ARRAY {1..64} OF CHAE; 

CBS : PACKED ARRAY (1..64) OF CHAR; 

H, C : INTEGER 

PROCEDURE BLANKLINE (N : INTEGER) ; 

VAR 

I ; INTEGER; 

EEGIN 

FOE I := 1 TC N DO WRITEIN (LISTING) ; 

END; 

PROCEDURE MAKEDOT; 

EEGIN 

WRITE (LISTING , ' *•); 

FOE C := 1 TO 128 DO BEGIN 

IF (C MOD 10) = 0 THEN WRITE (LISTING, '+’) ; 
ELSE WRITE (LISTING, ' ’) ; 

END; 

WRITE (LISTING, •*») ; 

WRITELN (LISTING) ; 

END; 

BEGIN (* CUT3AKMAP *) 
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CHS := ’01234567 89ABCDEFGHI JKLMNOPQRSTUVWXYZ " () <> 
(* keyboard dees not have all the characters 
that VAX-750 has *) ; 

FCR B := 1 TO 64 EO RCH {R} : = CHS{R}; 

E IA I KLINE (10) ; 

WRITELN (LISTING,’ REVISED GLOBAL REGION MAP'); 
EIA f KLINE ( 3) ; 

MAKEEOT; 

WRITE (LISTING,' ’); 

FOR C := 1 TO 64 £0 WRITE (LISTING, RCH {C} ) ; 
WRITELN (LISTING); 

MAKEEOT ; 

WITH BBEGM AP° DO EEGIN 
EIA l KLINE (3) ; 

MAKEEOT; 

FOR R := 1 TO 128 DO BEGIN; 

WRITE (LISTING, R: 3,'' ") ; 

FOE C := 1 TO 128 DO BEGIN 

WRITE (LI STING, RCH {BRMAP_PIX (R,C) MOD 64} ) ; 
ENE; 

WRITELN (LISTING); 

END; 

EIANKLINE (3); 

MAKEEOT; 

FOR R ;= 1 TO 128 DO BEGIN 
WRITE (LISTING, E:3, ' *); 

FOE C := 129 TO 256 DO EEGIN 

WRITE (LISTING, RCH {BRM AP_PIX {R, C} MOD 64} ) ; 
ENE; 

WRITELN (LISTING) ; 

END; 

EIA N KLINE (3); 

MAKEEOT; 

FCE E := 129 TO 256 DO EEGIN 
WRITE (LISTING, E;3, ’ ’); 



64 



FOE C := 1 TO 128 DO BEGIN 

WHITE (LISTING,RCH£ERMAP_?IX (R,C) MOD 64} ) ; 
END; 

WRITELN (LISTING) ; 

END; 

ElAiKIINE (3); 

MAKEEOT; 

ECE E := 129 To 256 DO BEGIN 
WRITE (LISTING, E: 3, » •); 

POE C := 129 TC 256 DO BEGIN 

WRITE (LISTING, ECH (BEMAP_PIX (R, C) MOD 64)) 
END; 

WHITELN (LISTING); 

END; 

MAKEDCT; 

END; 
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APPENDIX B 



PROCEDURE GEOMETRY 



PROCEDURE GEOMETRY 

(* This procedure calculates the area, 

perimeter, centroid, size and compactness 
for each of the region symbols used. 

Output format is for printer output. *) 

TYPE 

XYZ = ARRAY {1 . . 256, 1. . 256} OF INTEGER; 

ZYX = ARRAY { 1 . . 64} OF INTEGER; 

YZX = ARRAY {1.. 64.1. .2} OF INTEGER; 

VAR 



IE : 


XYZ; 










AREAP : 


ZYX; 










CENTER : 


YZX; 










PERIM : 


ZYX; 










SIZE : 


ARRAY 


{1. 


.64} 


OF 


REAL; 


COMPACT : 


ARRAY 


{1- 


. 64} 


OF 


REAL ; 


RCH : 


ARRAY 


{1. 


. 64} 


OF 


CHAR; 


CHS : 


PACKED 


ARRAY 


£1.. 


64} OF CHAR 



J, K, 1, R, C : INTEGER; 

PI : REAL; 

CONST 

MR = 256; 

MC = 256; 

NR = 64 ; 

NC = 4; 

PROCEDURE AREA 1 (IP : XYZ ; ISX : I STEGER ; ISY : INTEGER ; 

IRN; INTEGER; JF:ZYX ; NR: INTEGER) ; FORTRAN; 
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PROCEDURE CGBV1 (IP : XYZ ; IS X : I NT EGER ; IS Y : I NTEGER ; 

IR N : INTEGER ;JF: YZX ; NR : INTEGER) ; FORTRAN 
PROCEDURE PRMT1 (IP : XYZ ; IS X ;I NTEGER ; IS Y ; INTEGER ; 

IRN:INTEGER; JF:ZYX ; NR: INTEGER; 

NC : INTEGER) ; FORTRAN; 



BEGIN 

FOR R := 1 TO 256 DO BEGIN 
FOF C := 1 TO 256 DO BEGIN 

IP {R/C} := BREGMAP°.BRMAP_PIX{R,C} ; 

END: 

END: 

PI := 3.1415926; 

CHS := ’ 01 234567 8SABCDEFGHIJKLMNOPQR ST UVKXYZ H () <> 

(* other symbols missing on this keyboard *) ; 
FOR J := 1 TO 64 DO RCH {J} ;= CHS(J}; 

FCR K := 1 TO NR 10 BEGIN 
AREA1 (IP,MR, MC,K,AREAP,NR) ; 

IF AREAP (K) > 0 THEN BEGIN 
CGRV1 (IP,MR,MC,K, CENTER, NR) ; 

PRMT1 (IP,MR,MC,K,PERIM,NR,NC) ; 

SIZE £K) := 2 * AREAP {K} / PERIM {K} ; 

CCMPACT £K} := 4 * PI AREAP (K) / PERIM {K} ** 2; 

WRITELN (LISTING) ; 

WRITELN (LISTING) ; 

WRITELN (LISTING) ; 

WRITE (LISTING, ’REGION »,K:4); 

WRITE (LISTING,*, SYMBOL »,RCH}K MOD 64} : 1) ; 
WRITELN (LISTING); 

WRITE (LISTING, * AREA ’, AREAP {K} :5) ; 

WRITE (LISTING,* CENTROID AT », CENTER [K , 1} : 3, 
CENTER {K, 2} : 5) ; 

WRITELN (LISTING) ; 

WRITE (LISTING, ' PERIMETER ' , PERIM £K} : 5) ; 

WRITE (LISTING,*, SIZE »,SIZE{K}); 
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WRITS (LISTING,*, COMEACTNESS ', COMPACT {K} ) ; 
WEITELN (LISTING) ; 

ENE; 

END; 

ENE; 
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